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We study the effects of frustration in an antiferromagnetic film of FCC lattice with Heisenberg 
spin model including an Ising-like anisotropy. Monte Carlo (MC) simulations have been used to 
study thermodynamic properties of the film. We show that the presence of the surface reduces the 
ground state (GS) degeneracy found in the bulk. The GS is shown to depend on the surface in-plane 
interaction Js with a critical value at which ordering of type I coexists with ordering of type II. Near 
this value a reentrant phase is found. Various physical quantities such as layer magnetizations and 
layer susceptibilities are shown and discussed. The nature of the phase transition is also studied 
by histogram technique. We have also used the Green's function (GF) method for the quantum 
counterpart model. The results at low-T show interesting effects of quantum fluctuations. Results 
obtained by the GF method at high T are compared to those of MC simulations. A good agreement 
is observed. 
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I. INTRODUCTION 

Effects of the frustration in spin systems have been 
extensively investigated during the last 30 years. Frus- 
trated spin systems are shown to have unusual properties 
such as large ground state (GS) degeneracy, additional 
GS symmetries, successive phase transitions with com- 
plicated nature. Frustrated systems still challenge the- 
oretical and experimental methods. For recent reviews, 
the reader is referred to Ref. [l|. 

On the other hand, during the same period physics 
of surfaces and objects of nanometric size have also at- 
tracted an immense interest. This is due to important 
applications in industryi^i^i^i^ In this field, results from 
laboratory research are often immediately used for in- 
dustrial applications, without waiting for a full theoret- 
ical understanding. An example is the so-called giant 
magneto-resistance (GMR) used in data storage devices, 
magnetic sensors, ..^hIj^ In parallel to these experimen- 
tal developments, much theoretical effort has also been 
devoted to the search of physical mechanisms lying be- 
hind new properties found in nanometric objects such as 
ultrathin films, ultrafine particles, quantum dots, spin- 
tronic devices etc. This effort aimed not only at provid- 
ing explanations for experimental observations but also 
at predicting new effects for future experiments. ^'''^^ 

The aim of this paper is to investigate the effect of the 
presence of a film surface in a system which is known 
to be very frustrated, namely the FCC antiferromagnet. 
The bulk properties of this material have been largely 



studied as we will show below. In this paper, we would 
like to see in particular how the frustration effects on the 
nature of the phase transition in 3D are modified in thin 
films and how the surface conditions affect the magnetic 
phase diagram. To carry out these purposes, we shall use 
Monte Carlo (MC) simulations and the Green's function 
(GF) method for qualitative comparison. 
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The paper is organized as follows. Section II is devoted 
to the description of the model. We recall there proper- 
ties of the 3D counterpart model in order to better ap- 
preciate properties of thin films obtained in this paper. 
A determination of its GS properties is also given. In sec- 
tion III, we show our results obtained by MC simulations 
as functions of temperature T. The surface exchange in- 
teraction Js is made to vary. A phase diagram in the 
space (T, Jg) is shown and discussed. In general, the 
surface transition is found to be distinct from the tran- 
sition of interior layers. An interesting reentrant region 
is observed in the phase diagram. We also show in this 
section the results on the critical exponents obtained by 
MC multihistogram technique. A detailed discussion on 
the nature of the phase transition is given. Section IV is 
devoted to a study of the quantum version of the same 
model by the use of the GF method. We find interest- 
ing effects of quantum fluctuations at low T. The phase 
diagram (T, J^) is established and compared to that ob- 
tained by MC simulations for the classical model. Con- 
cluding remarks are given in section V. 



II. MODEL AND CLASSICAL GROUND STATE 
ANALYSIS 

It is known that the antiferromagnetic (AF) interaction 
between nearest-neighbor (NN) spins on the FCC lattice 
causes a very strong frustration. This is due to the fact 
that the FCC lattice is composed of corner-sharing tetra- 
hedra each of which has four equilateral triangles. It is 
well-known^ that it is impossible to fully satisfy simulta- 
neously the three AF bond interactions on each triangle. 

The analytical determination of the GS of systems of 
classical spins with competing interactions is a fascinat- 
ing subject. For a recent review, the reader is referred to 
Rcf. Ij. For the bulk FCC antiferromagnet, the Heisen- 
berg spins on a tetrahedron form a configuration charac- 
terized by two arbitrary angleSf^^ The ground state (GS) 
degeneracy is therefore infinite. This is also found in fully 
frustrated simple cubic lattice with classical Heisenberg 
spins;^ the GS is also characterized by two random con- 
tinuous parameters. To give an idea about the GS of the 
bulk FCC antiferromagnet,— let us imagine two planes, 
xz and where "0 intersects with the xz plane along 
the z axis and makes an angle (j) with the x axis. Two 
of the four spins make an angle 6 in the xz plane sym- 
metric with respect to the z axis. The other two spins 
make also the same angle, symmetric with respect to the 
z axis, but in the plane 'tp. It has been showni^ that 
the two angles and (j) are arbitrary between and tt. 
Note that when 6 = the spin configuration is coUinear 
with two spins along the -\-z axis and the other two along 
the —z one. The phase transition of the bulk frustrated 
FCC Heisenberg antiferromagnet has been studied;J^- In 
particular, the transition is found to be of the first order 
as in the Ising case i^^i^^'^^ Other similar frustrated anti- 
ferromagnets such as the HCP antiferromagnet show the 
same behavior,iS 

Let us consider a film of FCC lattice structure with 
(001) surfaces. To avoid the absence of long-range or- 
der of isotropic non Ising spin model at finite tempera- 
ture (T) when the film thickness is very small, i.e. quasi 
two-dimensional system2£, we add in the Hamiltonian an 
Ising-like uniaxial anisotropy term. The Hamiltonian is 
given by 

H = j.,,s.-s,~^A(^n' (1) 

where Si is the Heisenberg spin at the lattice site i, 
indicates the sum over the NN spin pairs Si and S j . 

In the following, the interaction between two NN sur- 
face spins is denoted by Js, while all other interactions 
are supposed to be antiferromagnetic and all equal to 
J = — 1 for simplicity. 

We first determine the GS configuration by using the 
steepest descent method : starting from a random spin 
configuration, we calculate the magnetic local field at 
each site and align the spin of the site in its local field. 
In doing so for all spins and repeat until the convergence 



2 




-1 -0.8 -0.6 -0.4 -0.2 J -I -0.8 -0.6 -0.4 -0.2 J 



FIG. 1: A ground state configuration of single plaquette (a) 
SI is of sublattice 1, (b) cos 612, (c) cos ^23, (d) cos 6^34. 
cos 6ij is cosine of angle between the two spins of sublattices 
i and j. 




(a) (b) 



FIG. 2: The ground state spin configuration of the FCC cell 
at the film surface: a) ordering of type I for Js < —0.5 b) 
ordering of type II for Js > —0.5. 



is reached, we obtain easily the GS configuration without 
metastable states. The result is shown in Fig. [T] 

We observe that there is a critical value = —0.5. 
For Js < Jg, the spins in each yz plane are parallel while 
spins in adjacent yz planes are antiparallel (Fig. [2^). 
This ordering will be called hereafter "ordering of type 
I" : in the x direction the ferromagnetic planes are an- 
tiferromagnetically coupled as shown in this figure. Of 
course, there is a degenerate configuration where the fer- 
romagnetic planes are antiferromagnetically ordered in 
the y direction. Note that the surface layer has an AF 
ordering for both configurations. The degeneracy of type 
I is therefore 4 including the reversal of all spins. 

For Jg > Jg, the spins in each xy plane is ferromag- 
netic. The adjacent xy planes have an AF ordering in 
the z direction perpendicular to the film surface. This 
will be called hereafter " ordering of type 11" . Note that 
the surface layer is then ferromagnetic (Fig. [Jb)- The 
degeneracy of type II is 2 due to the reversal of all spins. 

Without using a general methodji^ii^ let us calculate 
analytically the GS configuration in a simple manner for 
the present model. 

Consider a tetrahedron with the spins numbered as in 
Fig. [21 "S*!, 5*2, 5*3 and ^4 are the spins in the surface 
FCC cell (first cell). The interaction between and S2 
is set to be equal to Js (— I < < 0) and all others 
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are taken to be equal to J (< 0), and all Di = D for 
simplicity. The Hamiltonian for the cell is written as 

Hp = -i{8J,(Si-S2) + 8J(S3-S4) 

+ 6 J [Si • S3 + Si ■ S4 + S2 • S3 + Ss • S4] 

+ 2D [{Slf + {S^,f + {Sir + }• (2) 

Let us decompose each spin into two components: an xy 
component, which is a vector, and a z component Si = 
(s|,S'f). The numerical results shown above indicate 
that the spins have only z component. Taking advantage 
of this, we suppose that the xy vector components of the 
spins are all equal to zero. The angles 6i of Si with the 
z axis are then 

9i = 0, 63 =n 
'^2 — Pi, Oi—132- 

The total energy of the cell ([2]), with Si = ^, can be 
rewritten as 

Hp = + — + I — - - — cos /3i 1 cos ^1 

1 3 J 

+ - ( J — I5COS/32) cos/32 — ^ cos /3i cos /32 . (3) 

By a variational method, the minimum of the cell energy 
corresponds to 

dH D 

— ^ = JjSin^i + — cos^isin/3i 

dPi 2 



— sm Pi + — cos P2 sm Pi 



dHp 

df32 



3J 



cos (3i 



J D 



cos (32 



= 0, (4) 
sin/32 = 0. (5) 



The solutions of Eq. ^ and Eq. ([5]) corresponding to 
the minimal energy are 



cos/3i = — cos/32 = —1 for \ Js\ > 0.5, 
cos/3i = — COS/32 = 1 for \Js\ < 0.5. 



(6) 



Note that these solutions do not depend on D. The GS 
energy per spin is 



Hp = -D + J + Js for|J,|>0.5, 
Hp = -D + 2J - J, for \ Js\ < 0.5. 



(7) 



We see that the solution @ agrees with the numerical 
result. 



III. MONTE CARLO RESULTS 

In this paragraph, we show the results obtained by MC 
simulations with the Hamiltonian ([T|). The spins are the 
classical Heisenberg model of magnitude 5=1. 




FIG. 3: Magnetizations and susceptibilities of sublattices 1 
and 3 first two cells vs temperature for Js — —1.0 with L = 24 
and D = 0.1. Lj denotes the sublattice magnetization of layer 
j- 



The film size is L x L x where is the num- 
ber of FCC cells along the z direction (film thickness). 
Note that each cell has two atomic planes. We use here 
L = 12, 18, 24, 30, 36 and N, = 4. Periodic boundary 
conditions are used in the xy planes. The equilibrating 
time is about 10^ MC steps per spin and the averaging 
time is 2 X 10^ MC steps per spin. | J| = 1 is taken as 
unit of energy in the following. 

Before showing the results let us adopt the following 
notations. The sublattice 1 of the first cell belongs to 
the surface layer, while the sublattice 3 of the first cell 
belongs to the second layer. The sublattices 1 and 3 of the 
second cell belong, respectively, to the third and fourth 
layers. In our simulations, we used four cells, Nz — 4, i.e. 
8 atomic layers. The symmetry of the two film surfaces 
imposes the equivalence of the first and fourth cells and 
that of the second and third cells. It suffices then to 
show results of the first two cells, i. e. four first layers. 
In addition, in each atomic layer the two sublattices are 
equivalent by symmetry. Therefore, we choose to show 
in the following the results of the sublattices 1 and 3 of 
the first two cells, i.e. results of the first four layers. 

Let us show in Fig. [3] the magnetizations and the sus- 
ceptibilities of sublattices 1 and 3 of the first two cells, 
in the case where Js — —1. 

It is interesting to note that the surface layer has 
largest magnetization followed by that of the second 
layer, while the third and fourth layers have smaller mag- 
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netizations. This is not the case for non frustrated films 
where the surface magnetization is always smaller than 
the interior ones because of the effects of low-lying en- 
ergy surface-localized magnon modes. ^-'^"^^ One explana- 
tion can be advanced: due to the lack of neighbors surface 
spins suffer fluctuations due to the frustration less than 
the interior spins so they maintain their ordering up to 
a higher temperature. Let us decrease the Js strength. 
The surface spins then have smaller local field, so thermal 
fluctuations will reduce their ordering to a lower temper- 
ature. Figures m and [5] show respectively the cases where 
Js = —0.8 and -0.5. Near Js = —0.8 the crossover takes 
place: the surface magnetization becomes smaller than 
the interior ones for Js > —0.8. Note that the mag- 
netizations of second, third and fourth layers undergo a 
discontinuity at the transition temperature for Js = —0.8 
and -0.5. This suggests that the phase transitions for in- 
terior layers are of first order as it has been found for 
bulk FCC antiferromagnct.^^ This should be checked in 
the future. 

For weak |Jsl, there is only one transition for all lay- 
ers. An example is shown in Fig. [S] for Js = —0.1. 
Note that the first-order character disappears as there is 
no discontinuity of layer magnetizations at the transition 
temperature. 




FIG. 4: Magnetizations and susceptibilities of sublattices 1 
and 3 of first two cells vs temperature for Js = —0.8 with 
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FIG. 5: Magnetizations and susceptibilities of first two cells 
vs temperature for Js = —0.5 with L — 24 and D = 0.1. Lj 
denotes the sublattice magnetization of layer j. The suscep- 
tibility of sublattice 1 of the first cell is divided by a factor 5 
for presentation convenience. 




L = 24andD = 0.1. 
of layer j. 



Lj denotes the sublattice magnetization 



(b) 
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In the region —0.5 < Js < —0.45, there is an interest- 
ing reentrant phenomenon. To facilitate the description 
of this phenomenon, let us show the phase diagram in the 
space ( Js, Tc) in Fig. [T] In the region -0.5 < Js < -0.45, 



FIG. 6: Magnetization and susceptibility of first two cells vs 
temperature for Js = —0.1 with L = 24 and D — 0.1. Lj 
denotes the sublattice magnetization of layer j. 
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FIG. 7: Critical temperature vs Js with L = 24 and D — 0.1. 
Lj denotes data points for the sublattice magnetization of 
layer j. I and II denote ordering of type I and II defined in in 
Fig. (2] III is paramagnetic phase. The discontinued vertical 
line is a first-order line. See text for comments. 



FIG. 8: The magnetization M and the staggered magnetiza- 
tion A'Ist of first layer versus Js are shown, at T = 0.15, with 
L — 24 and D — 0.1. I and II denote ordering of type I and 
II defined in Fig. [2] III is paramagnetic phase. See text for 
comments. 



the GS is of type II as seen above. According to the phase 
diagram, we see that when the temperature increases 
from zero, the system goes through the phase of type II, 
undergoes a transition to enter the phase of type I before 
making a second transition to the paramagnetic phase 
at high temperature. This kind of behavior is termed 
as reentrant phenomenon which has been found by exact 
solutions in a number of very frustrated systemsj^i^l For 
a complete review on these exactly solved systems, the 
reader is referred to the chapter by Diep and Giacomini^^ 
in Ref. [l]. We note here that the reentrance is often 
found near the frontier where two phases coexist in the 
GS.i This is the case at Js = J^ = -0.5. 

The discontinued vertical hne at Jg — —0.5 is a first- 
order line separating phases I and II. The coexistence of 
these two phases which do not have the same symmetry 
explains the first-order character of this hne. To show 
it explicitly, we have calculated at T — 0.15 the magne- 
tization M and the staggered magnetization Mgt of the 
first layer with varying Js across -0.5. From the GS con- 
figurations shown in Fig. [21 M should be zero in phase 
I and finite in phase II, and vice- versa for Mst- This is 
observed at T — 0.15 as shown in Fig. [S] The large 
discontinuity of M and Mgt at Js = —0.5 shows a very 
strong first-order character across the vertical line in Fig. 

El 

Let us discuss on finite-size effects in the transitions 
observed in Fig. [31 to Fig. [51 This is an important ques- 
tion because it is known that some apparent transitions 
are artifacts of small system sizes. To confirm the ob- 
served transitions, we have made a study of finite-size 
effects on the layer susceptibilities by using the accurate 
MC multi histogram technique.— 

At this point, let us recaU that bulk Ising frustrated 
systems, unlike unfrustrated counterparts, have different 
transition natures: the antiferromagnetic FCC and HCP 
Ising lattices have strong first-order transitiou j^^i^^i^^ 
while the stacked antiferromagnetic triangular lattice has 
a controversial nature (see references in Ref. Wu). The 
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FIG. 9: Susceptibilities of layer 1 (left) and 2 (right) are shown 
for various sizes L as a function of temperature for Js = —0.1 
and D = 0.1. 



model studied here is the frustrated FCC film where sur- 
face effects can modify the strong first-order observed in 
its bulk counterpart. 

Our results show that transitions at Js = —1 and 
Js = —0.1 are real second-order transitions obeying some 
scaling law. Figure [9l shows the size effects on the maxi- 
mum of the susceptibilities of the first and second layers 
for Js = —0.1, while Fig. [TUl shows that of the third and 
fourth layers. As seen, the maximum of the susceptibili- 
ties -)^™^'^ increases with increasing L. 

Using the scaling law ^'^^^ oc L''^'^, we plot Inx"^^^ 
versus InL in Fig. [TTl The ratio of the critical expo- 
nents j/v is obtained by the slope of the straight line 
connecting the data points of each layer. 

Within errors the third and fourth layers have the same 
value oi^ /v which is neither 2D nor 3D Ising universality 
classes, 1.75 and 2, respectively. The same for the the 
values of the first and second layers. The exponent v 
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FIG. 10: Susceptibilities of layer 3 (left) and 4 (right) are 
shown for various sizes L as a function of temperature for 
Js = -0.1 and D = 0.1. 
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FIG. 12: The maximum value of ((InM)') = (E) - 
{M E) I (M) versus L in the In — In scale for Js — —0.1, where 
M is the sublattice order parameter. The slope of each line 
gives 1/v. Lj denotes the sublattice magnetization of layer j. 
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FIG. 11: Maximum sublattice susceptibility i^™'*" versus L 
in the In — In scale, for Js = —0.1 and D — 0.1. Lj denotes 
the sublattice magnetization of layer j. The slopes of these 
lines give the ratios of exponents "f/v. 



sublattice magnetizations and their derivatives for the 
first two layers in the case of Jg = —1. We find i^i — 
0.794 ± 0.022, 1/2 = 0.834 ± 0.027, 71 = 1.524 ± 0.0.040, 
and 72 = 1.509 ± 0.022. 



can be obtained as follows. We calculate as a function 
of T the magnetization derivative with respect to P — 
(kBT)-^: Vi = ((InAf)') = (E)-{ME) / (M) where is 
the system energy and M the sublattice order parameter. 
We identify the maximum of Vi for each size L. From 
the finite-size scaling we know that V{^^^ is proportional 
to L^/".^ We plot in Fig. [H In V/"'''^ as a function of 
InL for Js = —0.1. The slope of each line gives l/v. For 
the case Js = -0.1, we obtain 1/ = 0.822 ±0.020, 0.795 ± 
0.020,0.790 ± 0.020,0.782 ± 0.020 for the first, second, 
third and fourth layers. These values are far from the 
2D value (1/ = 1). We deduce 7 = 1.510 ± 0.010, 1.442 ± 
0.015, 1.412 ± 0.025, 1.395 ± 0.025. The values of ly and 
7 are decreased when one goes from the surface to the 
interior of the film. 

We show in Fig. [T3] and Fig. [TJ] the maximum of 



Let us discuss on the values of the critical exponents 
obtained above. These values do not correspond nei- 
ther to 2D nor 3D Ising models (72D = 1.75, 1^20 = 1, 
73D = 1.241, i^zD = 0.63). There are multiple reasons for 
those deviations. Apart from numerical precisions and 
the modest sizes we used, there may be deep physical 
origins. 

A first question which naturally arises is the effect of 
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FIG. 13: Maximum sublattice susceptibility ^'^'^^ versus L in 
the In — In scale, for Js = —1 and D — 0.1. Lj denotes the 
sublattice magnetization of layer j. The slopes of these lines 
give the ratios of exponents 'y/iy. 
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FIG. 14: The maximum value of ((InAf)') = (E) - 
(ME) I (M) versus L in the In — In scale for J a — —1, where 
M is the sublattice order parameter. The slope of each line 
gives l/i/. LkSj denotes one sublattice magnetization of layer 
j- 



the frustration. The 3D version of this model, as said 
above, has a first-order transition, with a very strong 
character for the Ising casei^ii^ii^ and somewhat less 
strong for the continuous spin models.— It has been 
shown that at finite temperature, the phenomenon called 
"order by disorder" occurs leading to a reduction of de- 
generacy: only coUinear configurations survive by an en- 
tropy effect The infinite degeneracy is reduced to 
6, i. e. the number of ways to put two AF spin pairs 
on a tetrahedron. The model is equivalent to 6-state 
Potts model. The first-order transition observed in the 
3D case is in agreement with the Potts criterion accord- 
ing to which the transition in q-state Potts model is of 
first-order in 3D for q > 3. 

In the case of a film with finite thickness studied here, 
it appears that the first-order character is lost. 



A first possible cause is from the degeneracy. Accord- 
ing to the results shown in the previous section, the GS 
degeneracy is 2 or 4 depending to Jg. If we compare to 
the Potts criterion according to which the transition is 
of first-order in 2D only when g > 4, then the transition 
in thin films should be of second order. That is indeed 
what we observed. 

Another possible cause for the second-order transition 
observed here is from the role of the correlation in the 
film. For second-order transitions, some arguments, such 
as those from renormalization group, say that the corre- 
lation length in the direction perpendicular to the film 
is finite, hence it is irrelevant to the criticality, the film 
should have the 2D character. If a transition is of first 
order in 3D, i. e. the correlation length is finite at the 
transition temperature, then in thin films the thickness 
effect may be important: if the thickness is larger than 
the correlation length at the transition, than the first- 
order transition should remain. On the other hand, if 
the thickness is smaller than that correlation length, the 
spins then feel an " infinite" correlation length across the 
film thickness. As a consequence, two pictures can be 
thought of: i) the whole system may be correlated and 
the first-order character is to become a second-order one 
ii) the correlation length is longer but still finite, the 
transition remains of first order. 

At this point, we would like to emphasize that, in the 
case of simple surface conditions, i.e. no significant devia- 
tion of the surface parameters with respect to those of the 
bulk, the bulk behavior is observed when the thickness 
becomes larger than a few dozens of atomic layers i^i^^^ 
surface effects are insignificant on some thermodynamic 
properties such as the value of the critical temperature, 
the mean value of magnetization at a given T, ... It 
should be however stressed that the criticality is very dif- 
ferent. It depends on the correlation length compared to 
the thickness: for example, we have obtained in the case 
of simple cubic films with Ising model the critical expo- 
nents identical to those of 2D Ising universality class up 
to thickness of 9 layersi^ Due to the small thickness used 
here, we think that the 2D character should be assumed. 

Now for the anisotropy, remember that in the case 
studied here, we do not deal with the discrete Ising model 
but rather an Ising-like Heisenberg model. The deviation 
from the 2D values may then result in part from a com- 
plex coupling between the Ising-like symmetry and the 
continuous nature of the classical Heisenberg spins. This 
deviation may be important if the anisotropy constant D 
is small as in the case studied here. 

To conclude this paragraph, we believe, from physi- 
cal arguments given above, that the critical exponents 
obtained above which do not belong to any known uni- 
versality class may result from different physical mecha- 
nisms. This is a subject of future investigations. 
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IV. GREEN'S FUNCTION METHOD 

We can rewrite the full Hamiltonian ([T]) in the local 
framework as 



l){S+S+ + SrSr) 



1 



icos0,, +1) {S+Sr +S-S+) 



+ I sin 9,, {S+ + Sr) - i sin (5+ + S^) 



^ 2 

2 



) <i> 



(8) 



where cos {9ij) is the angle between two NN spins. 

To study properties of quantum spins over a large re- 
gion of temperatures, there are only a few methods which 
give relatively correct results. Among them, the GF 
method is known to recover the exact results at very low- 
T obtained from the spin-wave theory. In addition, it is 
better than the spin- wave theory at higher temperatures 
and can be used up to the transition temperature with 
of course less precision on the nature of the phase tran- 
sition. It should be emphasized that the GF method is 
much better than other methods such as mean-field the- 
ories in estimating the value of the critical temperature. 
We choose here this method to study quantum effects at 
low T and to obtain the phase diagram at high T. 

The GF method can be used for non coUinear spin 
configurations In the case studied here, one has a 
coUinear one because of the Ising-like anisotropy. In this 
case, we define two double-time GF by^— 

G.,(i,0 - (9) 
F,,it,t') - «5r(t);5+(i')». (10) 

The equations of motion for Gij{t, t') and Fij{t, t') are 
written by 



i-G^At^t') = {[St{t),Sj{t')\)5{t-t') 

- {{[H,St{t)]-,S-{t'))), (11) 
^jF^Ai^t') - {[S;{t)^Sj{t')])6{t-t') 

(([7^,5^(^)];^-(^'))), (12) 

We shall neglect higher-order correlations by using the 
Tyablikov decoupling scheme^ which is known to be 
valid for exchange terms Then, we introduce the fol- 
lowing Fourier transforms 



gr.^r.'{u^,'k^y)e'^---^^^-^'\ (13) 



1 

A 



2tt 



+00 



fn.n' •) "^xy) ^ 



;k^y(R,-Rj) 



(14) 



where uj is the spin- wave frequency, kj;j, denotes the wave- 
vector parallel to xy planes, is the position of the spin 
at the site i, n and n' are respectively the indices of the 
layers where the sites i and j belong to. The integral over 
kajy is performed in the first Brillouin zone whose surface 
is A in the xy reciprocal plane. 

The Fourier transforms of the retarded GF satisfy a set 
of equations rewritten under the following matrix form 



M(cj)g = u. 



(15) 



where M [lo) is a square matrix {2Nz x 2Nz), g and u 
are the column matrices which are defined as follows 



and 



M(w) 



/ 51,"' \ 

fl,n' 




V / 




/ At 


Bi 


-Bi 


A-, 







(16) 



C 



c 



-D- 



A + 
—Bn^ 



Bn, 
An J 



(17) 



where 



UJ ± 



rl 



-J„(5^)(Z7) (cos0„ + l) 



J„(5^)Zcos0„-2/„ 

2Jn.n+l ('S'^+l) COS 0,^^ Jj^j^ 
2J„,n+l (5^+1) COS 0,^'') 



2J„,„-i(5,:_i)cos0i"i_i 



n,n4-l 
(a) 



.i(5,:_i)cos0i';Li 



B„ 



1 



J„(^^)(Z7) (cos( 



1) 



- Jn.n-l (Sn) (cOS 9l^J^_^ ± 1 

+ J„,„_i(5^)(cos^?i'Li±l 



(18) 

(19) 
(20) 
(21) 

in which, Z = 4 is the number of in-plane NN, O^^l-^j.^ 
the angle between two NN spins of sublattice 1 and 3 
belonging to the layers n and nil (see Fig. [2|), dn^n±i 



Jn,n+1 {Sn) 



cosC+i±l 
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the angle between two NN spins of sublattice 1 and 4, On 
the angle between two in-plane NN spins in the layer n, 
and 



1 



4 cos 



2 



cos 



Here, for compactness we have used the following no- 
tations: 

i) Jn and Dn are the in-plane interactions. In the 
present model J„ is equal to Jg for the two surface layers 
and equal to J for the interior layers. All Dn are set to 
be D. 

ii) Jn,n±i are the interactions between a spin in the 
n-th layer and its neighbor in the {n ± l)-th layer. Of 
course, J„.„_i = if n = 1, Jn.n+i — if n — Nz. 

Solving det|]V[| = 0, we obtain the spin- wave spectrum 
uj of the present system. The solution for the GF gn^n is 
given by 



9n 



Mn 

IMI 



(22) 



with |M|jj is the determinant made by replacing the n-th 
column of |M| by u in (fT6|) . Writing now 



(23) 



|M| = n - : 



one sees that lu^ (k-^y) 



Nz, are poles of the 



GF 

9n,n- (kajy) Can bc obtained by solving |]V[| — 0. 
In this case, gn,n can be expressed as 



gn,n — ^ ^ 

i 

where /„ {tOi {\i.xy)) is 

fn {^i (kxy)) 



[uJ - LOi {^^xy)) 



|M|„ {lo., (k.,)) 



(24) 



(25) 



Next, using the spectral theorem which relates the cor- 
relation function {S~S^) to the GF^^ one has 



lim ^ 

£^0 A 



dk^ 



(uj - is) ) 



+ 00 

-oo 

duj 



— (gn,n' (t^ + is) 

Ztt 



(26) 



where e is an infinitesimal positive constant and P = 
l/fc^T, ks being the Boltzmann constant. 

Using the GF presented above, we can calculate self- 
consistently various physical quantities as functions of 
temperature T. We start the self-consistent calculation 
from T = with a small step for temperature: 5 x 10^'^ at 
low T and 10~^ near Tc (in units of J/ks)- The conver- 
gence precision has been fixed at the fourth figure of the 
values obtained for the layer magnetizations. We know 
from the previous section that the spin configuration is 



collinear, therefore in this section, we shall use a large 
value of Ising anisotropy D in order to get a rapid nu- 
merical convergence. For numerical calculation, we will 
use D = A and J = — 1 and a size of 80^ points in the 
first Brillouin zone. 

Figure [15] shows the sublattice magnetizations of the 
first four layers. As seen, the first-layer one is larger than 
the other three just as in the case of the classical spins 
shown in Fig. [31 This difference in sublattice magnetiza- 
tion between layers vanishes at Jg ~ — 0.8 as seen in Fig. 
[TCI Again here, one has a good agreement with the case 
of classical spins shown in Fig. [4| 
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FIG. 15: Layer magnetization of first four layers vs temper- 



ature for Js 



-1.0 and D = 4. L, denotes the sublattice 



magnetization of layer j. Note that except the first layer, all 
other layer magnetizations coincide in this figure scale. 




FIG. 16: Layer magnetizations of first four layers vs temper- 
ature for Js = —0.8 and D = 4. Lj denotes the sublattice 
magnetization of layer j. Note that all layer magnetizations 
coincide in this figure scale. 



For Js > —0.8, the sublattice magnetization of the first 
layer is larger at low T and higher at high T as seen in Fig. 
[T7]for Js = —0.5. This crossover of sublattice magneti- 
zations comes from the competition between quantum 
fluctuations and the strength of Js'. when \Js\ is small, 
quantum fluctuations of the surface layer are small yield- 
ing a small zero-point spin contraction for surface spins 
at T = 0. So, surface magnetization is higher than the 
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interior ones. At higher T, however, small \Js\ gives rise 
to a small local field for surface spins which in turn yields 
a smaller surface magnetization at high T. This crossover 
has been found earlier in antiferromagnetic superlattices 
and films ,38'3£ 




iii) there is a reentrance between J, = —0.4 and Jg = 
—0.5. This is very similar to the phase diagram of the 
classical spins obtained by MC simulations shown in Fig. 

El 



0.5 1 




FIG. 17: Layer magnetization of first four layers vs temper- 



ature for Jg 



-0.5 and D = 4. L, denotes the sublattice 



magnetization of layer j. Note that except the first layer, all 
other layer magnetizations coincide in this figure scale. See 
text for comments on the crossover of surface magnetization. 



FIG. 19: Phase diagram obtained by the Green's function 
method with 7? = 4. Lj denotes the transition temperature 



of the sublattice magnetization of layer j. 
ments. 



See text for com- 



For Js = —0.1, there is no more crossover at low T as 
seen in Fig. [TBI Moreover, there is only a single transition 
at Tc ~ 2.65 for both surface and interior layers. 
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FIG. 18: Layer magnetization of first four layers vs temper- 
ature for Js = —0.1 and D = 4. Lj denotes the sublattice 
magnetization of layer j. Only at low T the surface magneti- 
zation is distinct from the other ones. 

We summarize in Fig. [19] the phase diagram for the 
quantum spin case obtained with the GF method. The 
vertical discontinued line indicates the boundary between 
ordered phases of types 1 and II. Phase III is paramag- 
netic. Note the following interesting points: 

i) for Js < —0.4 there is a surface transition distinct 
from that of interior layers, 

ii) for Js < —0.8, surface transition occurs at a tem- 
perature higher than that of interior layers. 



V. CONCLUDING REMARKS 

We have studied in this paper the properties of a thin 
film made from a fully frustrated material, namely the 
FCC antiferromagnet. We have considered both classi- 
cal and quantum Heisenberg spin model with an Ising- 
like single-ion anisotropy. The classical case was treated 
by Monte Carlo simulation while the quantum case was 
studied by the Green's function method. Several im- 
portant results are found in this paper. We found that 
the presence of a surface reduces the GS degeneracy of 
the fully-frustrated FCC antiferromagnet and there ex- 
ists a critical value of the in-plane surface interaction 
J^ ~ —0.5 which separates the GS configuration of type 
1 from that of type II. We have studied the phase transi- 
tion of the system. The surface spin ordering is destroyed 
in general at a temperature different from that of the in- 
terior layers. We found that in a small region just above 
Jg there is a reentrant phase: with decreasing T the sys- 
tem first changes from the paramagnetic phase to the 
type II phase, and then enters at a lower temperature 
into the type I phase. The critical behaviors of surface 
and interior layers have been shown and discussed. 

We hope that these unusual surface properties will help 
experimentalists to analyze their data obtained for real 
systems where frustration plays an important role. 

One of us (VTN) thanks the "Asia Pacific Center for 
Theoretical Physics" (South Korea) for a financial post- 
doc support and hospitality during the period 2005-2006 
where part of this work was carried out. 
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